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The presence of neutrino oscillations [1-2] leads to the question of possible CP violation 
in a (fully-coupled) three-family neutrino sector [3]. We suggest here a probe of CP that 
uses an effect of the modification of neutrino evolution within matter' [4-6] . The early 
two-family work of Refs. 4 and 5 arrives at a clear analytic result for propagation of 
neutrinos within a medium of constant density. Our method is also clearly illustrated in 
the two-family case. After development of our method we illustrate the three-family case 
with some numerical calculations. 

Neutrino propagation in material differs from propagation in vacuum because, in 
contrast to and V^, any Vg component in the neutrino beam can scatter from electrons in 
the material through charged current interactions. Consider neutrino passage through a set 
of N constant density material layers that we label 1 , 2, . . . , A'^. These could be materials of 
finite width or could be a set of labels that represent an (arbitrarily good) approximation 
to a material of continuously varying density. Because of multi-family transitions, the 
order of the materials matters. With n neutrino families, the n x n amplitude A^^ '^ is 

yKt)=A'^-''(t)xiX.O). (1) 
where y/ are the flavor states and t is the total time (or interchangably the total distance). 
The multiple-layer amplitude is the ordered product of the amplitudes for passage 
through the respective single layers, 

a12-^ = aUI..A^ (2) 

While Eq. (2) is general, to proceed further we want to start with two families, for 
which we know the single-layer amplitude is symmetric (see below). Then the 
amplitude for passage through the "reversed" material^, in which the order of layers is 
inverted, is 

A~ =A''-'' =A'' --A'A' = (a'A' ■■■A'^f 

(3) 

To see the consequences of this, recall that for two families with mass eigenvalues 
m,^ and a single (positive) mass difference factor dm^ = m2 - nii', and a Cabibbo matrix 
describing the flavor states in terms of the energy eigenstates in the usual way, the single- 
layer amplitude [4, 5] takes the general form 

7 1 

* * 

(In fact 77 is pure imaginary, but that is unimportant here.) We can then show by 
recurrence that the N layer amplitude has the structure 

a (3 

-p* a 

This follows because the amplitude forN+ 1 layers is 



(4) 



= A' 



(7) 



' In other work [3], the effects of matter imitate CP violation and it is necessary to distinguish the effects of 
matter from those of the presence of a CP-violating phase. 

^ We generally refer to quantities for the order 12. . .N as the "direct" case, and quantities for the order 
N...21 as the "reverse" case. 
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and this indeed has the structure of (7). It remains only to show by direct calculation 
starting from (4) that a two layer amplitude has the structure of (7), and this is easily 
done. This result can alternatively be derived in a continuum limit, for which the 
amplitude has an ordered exponential form. 

The combination of Eqs. (3) and (7) then tells us that 

'a -p* 

jS a 

Comparison between Eqs. (7) and (8) reveals that the probabilities Py = lAyl satisfy 

nreverse ndirect 

although only the diagonal elements of the amplitudes are equal. 

The presence of three families complicates the situation. Even if we could 
establish that Eq. (3) were valid, it would give only a useful relation for the diagonal 
matrix elements of the direct and reverse amplitudes; assuming the validity of (3) we find 

With this statement alone and unitarity, one can deduce the fact that for two families the 
off -diagonal probabilites are the same; for example. 
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For three families the off-diagonal probabilities for the direct and reverse probabilities 
would not generally be equal. Rather, for example, 
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Let us turn now to the question of three families, including the possibility that the 
mixing matrix is complex (CP in the mixing matrix). We begin by recalling the amplitude 
in vacuo. With the CKM matrix V connecting the flavor eigenstate l/^and energy 
eigenstate ^ according to the 1/^= VB„ one can easily show that 















A(r) = V 

























(Given the freedom to remove an overall phase, it is possible to use, for example, Ei = E 
+ nii l{2E), and pull out factors so that the central matrix gives only dependence on the 
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two mass difference factors dm2i = mi - mi and dm^i =m^ - mi of the problem.) 
Only if there is no CP violation, so that Vis real, is A sjmimetric. When matter is 
involved, and ignoring terms proportional to the unit matrix, the problem comes down to 
diagonalizing the matrix 
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This means finding a unitary matrix U such that 




























(11) 



The time evolution equation is then solved in the diagonal basis, and after transformation 
back to the flavor basis we have 


















U=U^D{t)U 



(12) 



(We remark here that if the consequences of a nonzero 5 are not to be lost, U cannot be 
found by using approximations based on a presumed neutrino mass hierarchy. Doing so 
leads to a decoupling in which the oscillation proceeds essentially through two-family 
steps, and this removes any effect of CP violation [7]. Similarly, if any dtrii^ = 0, then d 
can simply be set to 0.) 

At this point we can understand the effect of a CP-violating phase d in V. If 5=0, 
then y is real, and hence so is U, and A is symmetric. Equation (3) holds for the three 
channel case, and the respective diagonal elements of the direct and reverse processes are 
equal to one another. If, however, 5 is non-zero, then Vis no longer real, and neither is U. 
More specifically, we can generally write our primitive amplitude (12) with an index j for 
the material (including its width — t — > tj). Then for passage through a sequence of layers 
we have 

= ([// )* d; {-t,)u: (ui j d; {-t,)u; -(ui )* d; {-t.p; 

In other words, the more general form of Eq. (3) is 

^re.erse ^ ^j^^direct _^ ^ ^^direct _^ ^^^^ 

In particular, the diagonal matrix elements of the direct and reverse processes are no 
longer equal, and this provides a way to test for the presence of a nonzero value of 5. 

Of the diagonal processes available to us, we cannot use the process Vg Vg. That 
is because, given the fact that the charged current occurs in the electron sector, there is no 
CP violation visible in Vg — > Vg to leading order in the weak interaction [8]. We have 
verified the argument with, among others, direct numerical calculation with the 
parameters described below: the process Vg — > Ve remains independent of material order 
and indeed independent of 5. 

The process — > does allow us to test the effect described here. To illustrate 
we use a set of numbers that correspond to the so-called large-angle solar MSW scenario 
[9], namely ^32^ = ^31^ = 3.5 x 10"^ eV^ 5m2i^ = 5 x 10"^ eV^ sin^n = 0.10, sin023 = 
0.71, sin 012 = 0.53, with two layers, layer 1 corresponding to a material with density 8 
gm/cm^ (Vcc = 2.4 x 10"^^ eV) and layer 2 being that of the vacuum (or, for our purposes. 
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air). We consider three energies, £1 = 100 MeV, E2 = 500 MeV (corresponding to 
5ni'l{2E) = Vcc in material 1 for dni' = lO"'^ eV^), and £3 = 20 GeV (corresponding to a 
possible neutrino factory energy [3]). In Figs, la-c we plot for each energy the difference 
in probabilities for — > for the direct and reverse process. In the range from 5 = to 
5=1 rad, the 5 dependence is essentially linear in the magnitude of the differences, with 
no effect for 5 = to the largest effect for the largest 5. Aside from the height 
magnification there is no additional (5-dependence, and we restrict the plots presented 
here io 5= 0.5. The horizontal axis is the thickness of the air layer in units of 10^^ eV~^; 
we have in Figs, la-c taken equal amounts of vacuum and material. 

A well-known characteristic of the types of effects we are talking about is that 
they are compressed at low energies and stretched out at higher energies, and this is 
clearly visible in the three figures. We note the following: At £1, the effect tends to be 
more symmetric about the horizontal axis. Because in any hypothetical experiment the 
fine oscillations tend to be averaged over, the symmetry is a disadvantage. If one spreads 
the horizontal axis over a much larger range, say from to 100, then the envelope of the 
fine oscillations, which occurs in a series of "lobes" as is most visible in Fig. lb, moves 
slowly to either side of the horizontal axis, but clearly Ei is not the most interesting case. 
That is reserved for E2, where the lobes are clearly of one sign or another and the effect is 
reasonably large. The case Ej, also shows an average that is not symmetric about the 
horizontal axis, but over a larger distance scale. 

In Figure 2 we again consider the difference in probabilities for Vu for the 
direct and reverse process, but this time there is ten times as much air as material (Fig. 
2a) or ten times as much material as air (Fig. 2b); again the horizontal scale is the amount 
of air in the same units as in Fig. 1. The effect remains large and systematic. 

Finally we remark on an off-diagonal process, Vu Ve for definiteness. These 
processes have order dependence even in the absence of a CP- violating phase d. Figures 
3a-b plot the probability difference for the direct and reverse processes for 5 = and 5 = 
1 respectively. These probabilities certainly differ from one another; however, the 
standard for the detection of non-zero d is different from the diagonal case. 

Whether the phenomenon being described here is of eventual use will depend on 
as-yet unknown features of the neutrino spectrum. Detectable oscillations over 
managable distances would certainly render our results more relevant. It is nevertheless 
interesting that there is at least in principle a way to look for CP-violation in the neutrino 
sector without having to compare neutrino and antineutrino processes. 
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Figure Captions 

Figure 1 . Probability differences for finding v), as a function of time from production as a 
pure between the case where the neutrino beam passes through two different layers of 
equal thickness in one order and the case where the beam passes through the layers in the 
reverse order. The horizontal axis is the time for passage through (or thickness of) the air 
layer. Parameter values as described in the text, (a) E = 100 MeV; (b) E = 500 MeV; (c) 
E = 20 GeV. 

Fig 2 As in Fig. 1, except the thickness of the air and the material layer are different, (a) 
Air thickness 10 times that of material thickness; (b) air thickness 1/10 that of material 
thickness. 

Figure 3. Probability differences for finding Vg as a function of time from production as a 
pure v^i between the case where the neutrino beam passes through two different layers of 
equal thickness in one order and the case where the beam passes through the layers in the 
reverse order. The horizontal axis is the time for passage through (or thickness of) the air 
layer in units of 10^^ eV~\ Parameter values as described in the text, (a) 5=0; (h) 5= I. 
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